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Introduction 


A modem,  theoretical  investigation  of  the  nature  of  metallic 
cohesion  for  monovalent  atoms  necessitates  the  application  of  a form 
of  the  well  known  Schroedinger  equation  in  a determination  of  the  most 
stable  energy  position  for  the  valence  electron  while  considering  its 
relation  to  the  inner,  core  electrons  of  the  atom  as  well  as  to  the 
nucleus.  The  results  of  this  investigation  will  give  theoretical 
values  for  the  lattice  energy,  lattice  constant  and  compressibility 
of  the  metal  under  consideration.  This  method  was  applied  with  fair 
success  in  the  case  of  the  alkali  metals,  but  with  varying  success 
with  the  noble  metals. 

A comparison  of  the  results  of  work  done  by  K.  Fuchs  in 
connection  with  a theoretical  determination  of  the  lattice  energy, 
lattice  constant  and  compressibility  of  metallic  copper  , and  that 
done  by  R.  C.  Gunter  (unpublished)  with  metallic  silver,  brought  to 
light  the  fact  that  a fundamental  error  existed  in  one  or  the  other 
of  these  computations.  After  careful  consideration  the  writer  decided 
to  work  through  independently  the  calculations  for  both  of  these 
elements,  using  exiactly  similar  methods  in  each  case,  and  attempt  to 
determine  wherein  the  error  lay. 

The  final  results  indicated  the  error  to  be  with  the  previous 
calculations  for  copper,  the  writer  agreeing  fairly  well  with  Gunter's 
work  on  silver. 

The  writer  wishes  to  thank  Professor  Royal  M.  Frye  for  his 
unfailing  advice  and  encouragement  during  the  preparation  of  the  data 
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contained  in  this  dissertation,  and  Professor  Charles  0.  Ahonen  for 
his  suggestions  and  careful  reading  of  the  manuscript.  Finally,  he 
wishes  to  express,  however  inadequately,  his  appreciation  to  his  wife, 
Elizabeth  KacCullen  Merritt,  for  her  constant  help  in  the  transcribing 


of  numerical  data 
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Previous  Investigations  in  the  Field 


1.  E.  Wigner  and  F.  Seitz  on  Na 

Perhaps  the  first"*",  simplified,  quantum  mechanical  treatment  of 

cohesive  energy,  lattice  constant  and  compressibility  of  metallic 

2 

crystals  was  presented  by  E.  Wigner  and  F.  Seitz  in  connection  with 
the  constitution  of  metallic  sodium.  This  work  applied  the  free 
electron  assumption,  considering  that  the  electrons  in  the  K and  L 
shells  were  not  affected  by  a metallic  bond,  and  that  their  wave 
functions  were  the  same  as  in  the  free  state.  It  was  further  assumed 
that  the  wave  function  for  the  free  electron,  when  the  atom  was  in  a 
lattice  structure,  would  be  a repeating  function,  and,  because  of  the 
symmetry,  the  derivative  of  the  wave  function  at  every  cyrst alio graphic 
symmetry  plane  would  be  zero  perpendicular  to  this  plane.  This  vrauld 
give  the  initial  boundary  condition.  The  radial  Schroedinger 
(hereafter  abbreviated  S.)  equation  used  for  the  calculations  was  of 
the  following  form: 


It  should  be  pointed  out  here  that  this  is  actually  the  form  of 


'Previous  mathematical,  treatments  were  given  by  the  following: 

L.  Brillouin,  Die  Quant enstatistik,  Berlin  (1931) 

J.  C.  Slater,  Fhys.  Rev.  35,  509  (1930) 

E.  A.  Hylleraas,  Zeits.  f.  Physik  63,  771  (1930) 

H.  S.  Taylor,  H.  Eyring,  A.  Sherman,  J.  Chem.  Phys.  1,  68  (1933) 

"E.  Wigner  and  F.  Seitz,  Phys.  Rev.  43,  804  (1933),  46,  509  (1934) 

E.  Vagner,  Phys.  Rev.  46,  1002  (19347 
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the  Hartree  equation  for  a valence  electron.  The  potential  function 
should  include,  implicitly,  the  exchange  interaction  between  the 
valence  electron  and  the  core  electrons. 

The  potentials  used  by  Wigner  and  Seitz  were  obtained  from  the 
work  of  Prokofjew  , and  the  equation  was  solved  by  numerical  methods 
so  that  radii  corresponding  to  different  E values  could  be  determined. 
The  boundary  condition  could  be  applied  readily  so  that 

a _ d % _ 
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This  condition  would  occur  on  the  R-curve  at  points  where  a tangent 
would  cut  through  the  origin.  The  point  of  tangency  would  then  give 
the  r value  for  the  most  stable  energy  point.  In  general,  every  R 
curve  for  E less  (numerically)  than  that  for  the  most  stable  energy 
position  would  give  two  points  of  tangency.  For  E greater  (numerically) 
than  that  for  the  most  stable  energy  position  there  would  be  no  point 
of  tangency.  Curve  I shown  in  Fig.  1 (a)  shows  the  energy  curve 
obtained  by  this  procedure.  The  Fermi  energy  was  added  to  this,  giving 
curve  II  from  which  the  data  in  the  first  row  of  Table  I were  obtained. 

These  results  were  later  modified  by  a consideration  of  correlation 
energies  between  electrons  with  antiparallel  spin  so  that  the  data 
became  that  as  shown  in  the  second  row  of  Table  I.  In  this  instance 
it  was  believed  that  the  discrepancy  between  the  calculated  and 


Prokofjew,  Zeits.  f . P’nysik  58,  255  (1929) 


Rydberg  Units 


loO  2.0  3.0  4.0  5.0  6.0  7.0 

r,  Atomic  Units 


Table  I 


Theoretical  Data  for  Na 


1st  Determ.  25.6  Kg  Cal/mole 


Cohesive  Energy- 


Lattice  Constant  Compressibility 

4.20  A 1.6  x lO-11  cmVdyn 


2nd  Determ.  23.2 


4.75 


1.0  x 10"11 


Experimental  26.9 


4.23 


1.67  x 10-11 


experimental  values  could  be  accounted  for  by  a criticism  of  the 
expression  used  for  the  potential  within  the  s-sphere  at  the 
distance  r from  the  ion.  The  irork  in  this  second  determination 
was  actually  an  approximational  method  for  the  solution  of  a Fock 
equation  by  considering  exchange  interaction.  A discussion  of  the 
relation  between  the  Rartree  and  Fock  equations  will  be  given 
subsequently . 

2.  F.  Seitz  on  Li 

Following  the  method  adopted  in  the  work  on  the  constitution 
of  metallic  sodium,  Seitz^  undertook  the  solution  of  a similar 
problem  for  metallic  lithium.  A solution  of  Fock’s  equation  for 
the  free  electron  picture  gave  satisfactory  results,  but  correlation 
energies  had  to  be  obtained,  since  the  exchange  integral  in  the  Fock 
approximation  considers  only  electrons  of  parallel  spin.  In  this 
case  the  correlation  energy  was  borrowed  from  signer  and  Seitz' 
work  on  sodium,  and  was  employed  directly,  bringing  about  the  final 
results  shown  in  Table  II,  and  the  corresponding  energy  curve  shown 
in  Fig.  1 (b). 


^F.  Seitz,  Phys.  Rev.  47,  400  (1935) 
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Table  II 


Theoretical  Data  for  Li 


Cohesive  Energy  Lattice  Constant 

Calculated  33.8  Kg  Cal/mole  3.53  A 

Experimental  38.9  3.46 

3.  E.  Gorin  on  K 

The  work  of  Wigner  and  Seitz  on  Na  and  Li  was  published  in 
1933-34  and  1935,  respectively.  In  1936  Gorin'*  attempted  the 
determination  of  a set  of  Prokof jew  potentials  for  potassium,  but 
without  success,  probably  because  of  the  large  amount  of  interaction 

between  valence  and  core  electrons.  Accordingly,  he  applied  a set 

/ 

of  Hartree  functions0  which  were  available  for  the  core  electrons, 
and  evaluated  the  exchange  integrals  directly.  He  obtained  a value 
of  14.3  kilogram  calories  per  mole  as  against  the  experimental  value 
of  22.6  kilogram  calories  per  mole.  The  energy  curve  obtained  is 
reproduced  in  Fig.  1 (c). 

4.  K.  Fuchs  on  Cu 

The  first  attempt  at  the  application  of  Hartree  functions  to 

the  solution  of  a Fock  equation  for  one  of  the  noble  metals  was  made 

7 

by  K.  Fuchs  . He  applied  a set  of  available  Hartree  functions  for 
8 

copper  , and  evaluated  the  exchange  integral.  After  applying  the 

^E.  Gorin,  Physik.  Z.  Sowj.  9,  328  (1936) 

^D.  R.  Hartree,  Proc.  Roy.  Soc.  143.  506  (1934) 

7K.  Fuchs,  Proc.  Roy.  Soc.  151,  585  (1935) 

8D.  R.  Hartree,  Proc.  Roy.  Soc.  141.  282  (1933);  143,  506  (1934) 
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Table  III  Theoretical  Data  for  Cu  (preliminary,  after  Fuchs) 


Cohesive  Energy  Lattice  Constant  Compressibility, 
Caluulated  41.4  Kg  Cal/mole  r.  4.1  A 0.26  x 10“^-  cm' 


0.26  x 10“^-  cm  /dyne 
0.07 


Experimental  71.0 


3.6 


Fermi  energy  he  reported  the  results  shown  in  Table  III,  and  depicted 
by  the  energy  curves  in  Fig.  2.  Curge  (a)  is  without  exchange,  curve 
(b)  with  exchange,  and  curve  (c)  with  exchange  and  with  Fermi  energy 
added. 

Fuchs  then  claims  to  obtain  a better  value  for  the  compressibility 
by  considering  the  ion-ion  repulsion.  An  estimation  of  this  repulsion 
was  made  from  the  following  expression: 


n = number  of  nearest  neighboring  atoms  (12  in  this  instance) 
/°  = inter-atomic  distance 


(r)  = effective  nuclear  charge  (Hartree's  Z.  ) 


A ty4rr  = negative  charge  density  of  an  ion  (A  is  the  Laplacian) 
r^  and  r2  are  distances  from  two  neighboring  nuclei 


kK  - 3/5  e1/3  a0  (3n/f*S)Z/3 


kA  = 3/2  e2/3  aQ  (3/32rr2)1/3 


aQ  = Bohr  atomic  unit  of  length 
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Table  IV  Ion-Ion  Repulsion  Energy  for  Cu 
p (Atomic  Units)  4.5  5.0  5.5 

Wg  (Rydberg  Units)  0.197  0.0957  0.0161 

The  results  are  shox^m  in  Table  IV.  When  this  energy  is  considered, 
according  to  Fuchs,  curve  (d)  of  Fig.  2 is  obtained,  and  the  data 
are  as  shown  in  Table  V. 

Table  V Theoretical  Data  for  Cu  (final,  after  Fuchs) 

Cohesive  Energy  Lattice  Constant  Compressibility 
Calculated  34.1  Kg  Cal/mole  4.2  A 0.069  x 10”-*--1-  cm2 /dyne 

Experimental  71.0  3.6  0.070 

A criticism  of  these  results  together  with  a completely  new  set 
of  calculations  vail  form  a major  portion  of  this  dissertation. 

5.  R.  C.  Gunter  on  Ag 

It  appears  that  nothing  more  was  done  in  this  field  following 

9 

Fuchs'  work  until  Gunter  attacked  a similar  problem  for  silver. 

He  applied  exactly  the  same  method  that  Fuchs  did  for  copper,  a set 
of  Hartree  functions’^  having  been  available.  The  Fermi  energy  was 
added  to  the  energy  curve  obtained  by  solving  a Fock  equation  with 
the  exchange  integral.  The  final  results  are  shown  in  the  data  of 
Table  VI  and  the  curves  of  Fig.  3.  These  curves,  (a)  and  (b), 

q 

7R.  C.  Gunter,  Ph.D.  Dissertation,  Boston  University  (1942) 

1(k  M.  Black,  Kern,  and  Proc.  Manchester  Lit.  Phil.  Soc.,  59,  No.  3, 

28  (1934-35) 
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represent,  respectively,  energy  obtained  -without  and  with  the  Fermi 
energy  added.  No  ion-ion  repulsion  or  anti-parallel  spin  correlation 
was  considered. 


Table  VI 
Calculated 


Theoretical  Data  for  Ag  (after  Gunter) 

•Cohesive  Energy  Lattice  Constant 

$5.3  Kg  Cal/ mole  2.72  A 


Experimental 


64.5 


4.03 
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The  Mature  of  the  Problem 


Both  copper  and  silver  form  face  centered  crystal  lattices. 

With  similar  crystallographic  structures  these  metals  might  be 
expected  to  possess  relative  values  of  such  physical  quantities  as 
lattice  constant  and  compressibility  which  would  be  in  accordance 
with  the  relative  sizes  of  the  individual  atoms.  Actually,  experimental 
evidence  shows  both  the  lattice  constant  and  compressibility  of 
silver  to  be  greater  than  that  of  copper.  Therefore,  a consistent, 
theoretical  picture  of  metallic  copper  and  silver  should  show  the 
most  stable  energy  point  of  the  former  to  be  nearer  the  nucleus  than 
that  of  the  latter.  This  would  make  the  lattice  constant  and 
compressibility  of  silver  greater  than  that  of  copper,  consistent 
with  experiment. 

A careful  study  of  the  energy  curves  obtained  for  copper  by 
Fuchs  and  for  silver  by  Gunter  shows  the  stable  energy  point  to  be 
higher  for  copper  than  for  silver.  Obviously,  if  the  same  method 
had  been  used  on  both  atoms,  one  of  these  must  be  in  error.  Also, 

Fuchs  reports  that  the  reason  for  considering  the  ion-ion  repulsion 
was  to  obtain  a better  value  for  the  compressibility.  Actually,  the 
value  he  would  have  obtained  without  considering  the  ion-ion  repulsion 
would  have  been  better,  since  the  addition  of  the  ion-ion  repulsion 
energy  gives  a larger  value  for  rQ,  the  stable  energy  point,  and  hence 
a higher  value  for  the  compressibility. 

—1  1 2 

Fuchs  reports  the  compressibility  to  be  0.26  xlO  x cm  /dyne 

—1  "I  p 

without  ion-ion  repulsion,  and  0.069  xlO  cnr/dyne  with  ion-ion 
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repulsion  - - a definite  decrease.  In  this  connection,  if  the  curves 
he  reports,  and  which  were  reproduced  in  Seitz1  book,  The  Modern 
Theory  of  Solids,  were  used  he  could  not  have  possibly  obtained  the 
compressibility  values  he  claims  to  have  obtained.  Actually,  these 
values,  as  obtained  from  his  curves,  should  be  approximately 
0.32  x 10“^"  cm2/dyne  and  0.545  x 10“^  cm^/dyne  without  and  with 
ion-ion  repulsion,  respectively. 

In  view  of  this  discrepancy,  the  writer  decided  to  work  through 
independently  the  wave  functions  without  and  with  exchange  for  both 
copper  and  silver,  and  determine  if  possible  the  source  of  error. 

If  the  theoretical  stable  energy  position  for  copper  fell  closer  to 
the  nucleus  than  that  of  silver,  and  if  the  values  for  the  lattice 
constants  and  compressibilities  were  of  the  proper  order  of  magnitude, 
it  was  felt  that  the  investigation  would  have  been  well  worth  the 


effort 


Theoretical  Background  and  Approach 


1.  The  Schroedinger  Equation 

The  cellular  raethod  of  division  of  the  atoms  in  the  solid  state 
is  applied  to  the  face-centered  crystal  lattices  of  copper  and  silver. 
Each  atom  is  thus  enclosed  by  a polyhedron  which  can,  to  a first 
approximation,  be  considered  a sphere.  It  is  assumed  that  the  valence 
electrons  will  occupy  the  most  stable  energy  positions  between  atoms. 

Thus  the  -function  will  be  a repeating  function  of  the  form 

-w  cTT a K*  fT 


where  r is  the  position  vector,  and  y/r)  has  the  periodicity  of 
the  lattice,  k is  a wave  number  vector  arising  in  the  solution  of 
the  S.  equation,  and  having  the  components 

kx  = 1/Nxa  » = m/Nyb  » kz  = n/Nzc  » 

where  1,  m and  n are  arbitrary  integers,  N_,  and  Nz  are  the  numbers 

of  cells  in  the  respective  directions,  and  a,  b,  and  c are  the  inter 

atomic  distances  in  the  respective  directions.  The  space  defined  in 

this  manner  is  commonly  referred  to  as  k-space.  The  -function 

will  repeat  from  atom  to  atom,  and  will  be  smooth  and  constant 

between  atoms  at  the  stable  energy  point. 

The  5.  equation  to  be  used  in  this  investigation  will  be  of  the 
form  for  a single  electron,  and  the  potential  and  exchan0e  functions 
will  be  inserted  in  numerical  form  in  the  process  of  solution  of  the 
equation. 

We  shall  start  with  the  equation  in  the  radial  form 


« 

. 


. 


14 


dr* 


Srr^m.  . 

+ — r-(E  - V)  R = 0 


where  R = r 4>  (r),  r being  the  distance  out  from  the  center  of  the 
atom,  and  being  the  eigenfunction  for  the  valence  electron.  It 
may  be  well  to  point  out  here  that  this  equation  is  also  a reduced 
form  of  the  Kartree  equation  for  a single  valence  electron.  No 
exchange  between  valence  and  core  electrons  is  considered  in  the 
solution  of  this  equation  as  it  stands. 

If  we  use  atomic  units  for  length  and  energy: 

1 atomic  unit  of  length  = h^/4rrmiie^  = 0.528  A 
1 atomic  unit  of  energy  = 2n^meVh^  = 1 Rydberg  unit 
1 Rydberg  unit  = 13*53  electron  volts/ atom 
= 312  kilogram  calories/mole 
If  we  now  set  E = - € /2  our  equation  becomes 

^ + (2V  - € ) R = 0 


which  is  the  form  in  which  it  is  used  for  determining  the  R-function 
without  exchange. 

2.  The  Exchange  Integral 

The  energy  involved  in  the  interaction  between  the  valence 
electron  and  the  inner  core  electrons  is  contained  in  what  is  called 
the  exchange  integral.  This  was  used  by  Fuchs  in  his  calculations 
for  copper,  and  occurs  as  an  added  term  to  the  potential  energy  function 
in  an  S.  equation  to  give  a Fock  equation  for  a valence  electron 
which  takes  into  account  exchange  interaction  for  parallel  spins 
explicitly  so  that  this  is  not  required  of  the  coulomb  potential 


. 


, 


. 

- 
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function.  The  general  form  of  the  exchange  integral  is 


J go,aO  RgvmjV 


where 


GC*.jV)  = ^TV^-rx 


'«  J i„-,i  - 


*Y  being  the  angle  between  r and  r1 . The  summation  is  over  all  the 
inner  electrons  with  the  same  spin  as  the  valence  electron  so  would 
include  half  the  number  of  electrons  in  each  sub-shell.  The  's 
are  of  the  form  » S.  r \ 

_J Ta  U)  . 

^ \Qy&,V) 


%^zf*  X> 

being  a spherical  harmonic  normalized  by  the  following 

condition 

J7Y/Y,  ^ 


Thus  the  valence  electron  function  for  the  lowest  state  is 

U) ! — t (a) 

t ~ ajtt  s\ 

Following  Fuchs,  G(r, r*)  may  be  put  into  the  form 


where 
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11  r'  r 


if  r«  > 


The  complete  Fock  form  of  the  equation  for  the  valence  electron  would 

+ 5^=2-(e-  v(a))R  RoOJLa 

<A  J\X  o 


then  be 


. 


. ' 


1 


1 , < 

1 


1 


By  employing  atomic  units,  as  previously,  this  equation  reduces  to 

-v-  (tv-  e)  ^ -v  2 f G-Cn.jv)  Roo  i-*'  = o 

Ir*2-  o 

A complete  derivation  of  the  Hartree  and  Fock  equations  is  given  by 

Seitz'*’"'’  and  -will  not  be  dealt  with  here.  However,  since  the  Hartree 

12 

functions  as  obtained  by  the  method  of  the  self  consistent  field 
are  employed  in  the  calculations  for  both  copper  and  silver,  a brief 
mention  of  this  method  is  not  out  of  place. 

It  is  assumed  that  the  electronic  Vp  -function  for  the  entire 
atom  may  be  represented  as  the  product  of  the  individual  electron 
-functions: 

^ <*.  M’jd.-wO 't' 


Now  each  should  satisfy  an  S.  equation  of  the  form 

3 


A 


^ IV 


where  the  bracket  expression  for  the  total  potential  takes  into 
account  the  effect  of  the  neighboring  electrons  as  well  as  the 
nucleus,  and  6^  is  the  energy  parameter.  The  summation  is  over  all 
i except  i = j. 

For  the  simple  non-trivial  case  of  helium,  the  S.  equations  for 
the  individual  electrons  become 

-V-  = ° 


and 


-v  25^2. [e  - vtcm]M>r  =° 
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F.  Seitz,  Modern  Theory  of  Solids,  677-680,  McGraw-Hill  (1940) 
D.  R.  Hartree,  Proc.  Cambridge  Phil.  Soc.  24,  105  (1928) 
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- « 


ft 


Here  V-^  and  are  the  total  potentials.  Now,  V-j  is  estimated,  and  the 
first  equation  is  solved  for  . Then,  since 

V ;oo  ^ + czf  1^1  A 

SX  J 

Vn  may  be  determined. 

Using  this  value  for  V^,  the  second  equation  may  now  be  used  to 
solve  for  . From  this,  using  the  expression 

J rv»  , z 


V,(n)  = + * 


we  may  obtain  a new  approximation  to  V^.  This  is  continued  until 


V* 


V, 


- > 


V,, 


is  less  than  any  desired  ratio,  at  which  point 


the  assumed  and  calculated  fields  are  said  to  be  consistent  with 
each  other.  The  V is  usually  reported  in  -terms  of  effective  nuclear 
charge,  Z^  = rV,  and  the  difference  between  any  two  Z^  may  be  made  as 
small  as  one  wishes.  The  effective  nuclear  charge  is  that  charge 
which  must  be  placed  at  the  nucleus  position  to  give  a field  at  the 
distance  r which  is  the  same  as  the  calculated  field. 


3.  The  Fermi  Energy 

The  mean  kinetic  energy,  called  the  Fermi  energy,  of  a gas 
obeying  Fermi-Qirac  statistics  is  given  in  Rydberg  units  by 

aJlf'  l . 'z./clca 


W,  = , 

T \ O 


( 3 V/3  \ 
UtrJ  \, 


V V5 


where  V is  the  atomic  volume  equal  to  4/3  Trr^.  Fuchs  showed  that  if 


the  electrons  are  not  considered  "free"  the  more  general  expression 

- 0.10015 

1*  o 


should  be 


Vir.-~v'v<..  A-v( 


where 


g being  the  vector  drawn  from  the  origin  in  k-space  to  the  nearest 
point  of  the  first  zone  boundary  for  the  face  centered  cube.  He 
further  showed  that  this  expression  gave  results  very  close  to  that 
given  by  a free  electron  distribution.  Consequently,  for  all  Fermi 
energy  calculations  in  this  work  the  simpler  formula  will  be  used. 

4.  The  Cohesive  Energy 

If  EQ  is  the  self  energy,  that  is,  the  energy  without  exchange, 
of  the  valence  electron,  E^  the  Fermi  energy,  and  E.  the  negative 
of  the  ionization  energy,  we  can  express  the  cohesive  energy  as 


The  cohesive  energy  is,  then,  the  difference  between  the  ionization 
energy  and  the  energy  the  electron  has,  including  Fermi  energy,  when 
the  atom  is  in  a lattice  structure. 

5.  The  Compressibility 

This  quantity  may  be  obtained  theoretically  from  the  expression 


If  it  is  assumed  that  all  of  the  curvature  in  the  E vs.  r curve 
arises  from  the  Fermi  energy,  and  in  the  case  of  both  copper  and 
silver  this  is  true,  the  curve  with  exchange  but  without  Fermi  energy 
being  approximately  linear  in  this  region,  then  the  compressibility 
expression  reduces  to 


E.  = E + E„  + E. 
c o f 1 


160rrnm,  . j/s  ^ _n  ^ o 

K = ~9"2'"'(2rT/3)  / r0^  = 0.00190  x 10  J"L  cm2/dyne. 


This  is  a direct  fifth  power  relationship  between  the  compressibility 


' • ''it 
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and  the  stable  energy  point.  It  is  obvious  that  the  compressibility 
might  also  be  related  directly  to  the  lattice  constant  by  a similar 
fifth  power  expression. 

The  experimental  values  of  the  lattice  constant  and  compressibility 

13 

were  available  for  both  the  alkalis  and  noble  metals  . Because  of 
the  relationship  mentioned  above,  a log-log  plot  of  the  experimental 
values  of  compressibility  against  rQ  should  show  the  reliability  of 
the  theoretical  compressibility  expression.  Figure  4 shows  a plot 
of  compressibilities  against  experimental  rQ  for  both  the  alkalis  and 
noble  metals.  The  heavy  line  drawn  is  the  theoretical  curve,  ^t  can 
be  seen  that  potassium  is  the  only  one  of  these  elements  which  falls 
directly  on  this  curve.  For  the  others,  agreement  with  both  lattice 
constant  and  compressibility  values,  as  calculated  theoretically, 
would  be  impossible.  Gold  should,  theoretically,  fall  above  copper 
and  silver,  but  instead  is  lower  than  either  in  compressibility.  The 
reliability  of  the  fifth  power  expression  for  obtaining  theoretical 
values  of  the  compressibility  from  theoretical  values  of  rQ  might 
therefore  be  open  to  question. 

6.  The  Lattice  Constant 

The  lattice  constant  is  calculated  from  the  familiar  expression 

1/3 

for  a face  centered  crystal  , a = (l6rr/3)  (0.528)ro.  This  gives 

the  value  in  Angstrom  units  when  r is  in  atomic  units,  the  factor 
0.528  being  the  conversion  factor. 

13 
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The  Numerical  Method  Used  in  Attacking 
the  Problem 


This  problem  involves  the  solution  of  two  differential  equations 
and  one  integral  by  numerical  methods  for  each  of  the  two  metals 
under  consideration.  The  equation  without  exchange  is  of  the  form  of 
a simple  S.  equation, 

^2r 

jps  ♦ (2v  - e )R  - 0. 

The  equation  'with  exchange  is  of  the  following  form: 


+ (2V  - € )ft 
dr 


+ 2fe2Zl 

K 0 


V 0(r*  )R(r»  )dr' 


M 


'O 

• The  differentiation  involved  in  the  solution  of  the  first  of 

these  equations  was  accomplished  by  the  method  of  finite  differences, 

The  equation  was  first  put  in  the  following  form: 

2 

A^R 


2^2  + (2V  - €)R  = 0 


A table  was  then  constructed  similar  to  that  given  below: 

AR 


2 

A^R 


R, 


Ri 


R. 


S1  - ao 
*2  ' \ 


R2  ™ 2Rl  + Ro 


R. 


n-2 


R, 


n-1 


n 


Vi  - V2 

lRn  ~ Rn-1 


^ “2Rn-l  + Rn-2 


. 


J 


- 


Substituting  the  general  quantities  in  the  equation  we  obtain,  after 
rearranging, 

Rn  = 2V!  - V2  - 8n-l(^)(2V  - 6 ) 

Values  of  Rjj  may  be  determined  by  inserting  values  of  R^  n and  Rj^ 
together  with  V and  6 . The  function  is  known  to  go  to  zero 

at  the  origin,  and  any  initial  slope  may  be  taken  at  the  offset 
providing  the  £r  value  is  small  enough.  For  small  r values  £r  was 
taken  to  be  0.02  atomic  unit.  This  was  increased  as  the  R-curve 
became  less  fluctuating.  R-curves  for  different  energy  parameters, 

€ = - 2E,  were  obtained  by  solving  this  equation  for  each  € used. 

An  exactly  similar  method  was  employed  to  obtain  the  R-curve 
with  exchange,  the  equation  used  being 

Rn  - 2Rn-I  - V2  * Vl^^2v  ~ 6 ) - 2(Ar2)  J- 
where  ^ is  the  exchange  integral  which  is  a function  of  r. 

The  integration  involved  in  the  exchange  integral  was  performed 

by  first  obtaining  values  of  R(r')  ^ K f (r)  fKn(r’)  where  the 

K,*  * 

summation  is  over  ell  electrons  having  the  same  spin,  and,  as  mentioned 
previously, 

f -* 

I Pv 

if  r*  4 r 


fi-  I 



K 


rv 


SL 


if  r'  ^ r 


V JV 

The  f’s  are  the  radial  wave  functions  for  the  individual  electrons 

14  1 

as  worked  out  by  Black  for  silver  and  Hartree  J for  copper.  It 


14, 


See  reference  10 


15 


See  reference  8 


might  be  mentioned  that  these  tables  also  include  data  for  Z used  in 

Jr 

the  determination  of  the  values  for  the  coulomb  potential  V.  To 

illustrate  the  meaning  of  the  f's,  f^  q(0.02)  would  represent  the  value 

of  the  radial  wave  function  for  a 3s  electron  at  a distance  of  0.02 

atomic  unit  from  the  center  of  the  atom.  A complete  table  was  made 

for  both  silver  and  copper,  showing  the  results  of  the  summation  for 

all  possibilities  from  r = 0.02  to  r = 10  atomic  units.  Using  copper 

as  an  example,  the  actual  expansion  of  the  summations  to  be  made  for 

obtaining  the  value  of  the  exchange  integral  to  be  inserted  in  the 

differential  equation  for  a given  r is  shown  in  Table  VII  for  r = 0.20 

atomic  unit.  A summation  must  be  performed  for  each  r*.  The  relation 

of  each  term  of  the  summation  to  the  various  electronic  energy  states 

is  indicated  by  the  notation  Is,  2s,  2p,  etc.  at  the  top  of  the  table. 

The  expression  ^K^f(r)f(r')  is  called  A'r  r,  and  is  normalized 

by  the  factor  1/Arr  so  that  values  were  obtained  for  A*  ,^,/i+rr. 

r,  r 

Complete  sets  of  A'^  ^,/Att  values  are  shown  in  Tables  VIII  and  IX 
for  Cu  and  Ag,  respectively.  Since  these  tables  are  symmetric  about 
the  diagonal,  only  those  values  lying  on  and  above  the  diagonal  are 
recorded.  These  values  must  be  multiplied  by  the  corresponding 
R-functions  obtained  by  solving  the  S.  equation  without  the  exchange 
integral.  The  results  of  these  products  are  plotted  in  Fig.  5. 

This  is  a plot  of  A_  . against  r'  where  A , = Rfr')!'  . = 

r,r'  » r,r»  r,r' 

H(r')  > K.f.  .(r)  f„  .(r1).  The  area  under  this  curve  is  determined 

K,J>  * K>* 

mechanically,  and  represents  the  value  of  the  exchange  integral  for 


the  given  r. 
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.70  5(-2.15)(-.1070)  4.29(2.04)0.0920  5(-. 638)1. 271  4.29(0.430) (-1.206)  2.04(0.349)0.920 

.80  5 (-2. 15) (-.0474)  3.75(2.04)0.0425  5(-. 638)1.197  3.75(0.430)(-1.118)  1.562(0.349)0.890 

1.0  5 (-2. 15) (-.00837)  3(2.04)0.00836  5(-.638)0.835  3(0.430)(-.863)  1(0.349)0.797 
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Table  VIII 

Values 

of  A'/4rr  for 

Cu 

r .02 

.04 

.06 

.08 

.10 

.12 

r' 

.02  4,770 

5,360 

4, 530 

3,430 

2,430 

1,710 

.04  , 

3,010 

2,550 

1,925 

1,368 

932 

.06 

1,446 

1,093 

776 

531 

.08 

624 

),ij, 

-r-r-r 

274 

.10 

260 

178.2 

.12 

106.5 

r .14 

.16 

.18 

.20 

.30 

.40 

r 1 

.02  1,098 

712 

457 

287 

18.93 

-2.73 

.04  621 

403 

259 

163.6 

11.89 

-1.108 

.06  354 

238 

149.8 

95.5 

8.32 

-0.0589 

.08  205 

135.2 

88.8 

57.8 

6.56 

0.622 

.10  121.2 

81.3 

54.6 

36.5 

5.64 

1.045 

.12  73.8 

50.8 

35.2 

24.6 

5.12 

1.297 

.14  55.9 

33.4 

24.2 

17.66 

4.77 

1.440 

.16 

23.3 

17.66 

13.56 

4.49 

1.517 

.18 

13.69 

11.01 

4.33 

1.549 

.20 

9.28 

3.97 

1.557 

.30 

2.82 

1.549 

1.639 
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Table  VIII  (Cont.)  Values  of  k'/l&  for  Cu 


r .50 

.60 

.70 

.80 

1.0 

1.2 

r 1 

.02  -0.270 

0.990 

1.466 

1.591 

1.210 

0.817 

.04  -0.1356 

0.385 

0.594 

0.657 

0.507 

0.344 

.06  -0.0599 

-6.67xlO-3 

0.0356 

0.063 8 

0.0609 

0.0374 

.08  -7. 96x10" 5 

-0.239 

-0.301 

-0.294 

-0.208 

-0.1314 

.10  0.0536 

-0.360 

-0.481 

-0.486 

-0.353 

-0.230 

.12  0.1115 

-0.393 

-0.558 

-0.567 

-0.415 

-0.271 

.14  0.1750 

-0.371 

-0.545 

-0.570 

-0.420 

-0.276 

.16  0.246 

-0.305 

-0.489 

-0.522 

-0.390 

-0.256 

.18  0.324 

-0.213 

-0.366 

-0.424 

-0.336 

-0.221 

.20  0.408 

-0.1039 

-0.292 

-0.342 

-0.267 

-0.1774 

.30  0.363 

0.505 

0.317 

0.212 

0.1097 

0.0772 

.40  1.263 

0.989 

0.777 

0.619 

0.393 

0.235 

.50  1.534 

1.278 

1.035 

0.839 

0.519 

0.321 

.60 

1.413 

1.146 

0.924 

0.572 

0.352 

.70 

1.165 

0.933 

0.574 

0.352 

.80 

0.909 

0.555 

0.340 

1.0 

0.486 

0.297 

1.2  0.258 


Table  VIII  (Cont.)  Values  of  A */4n  for  Cu 


r 1.4 

1.6 

1.8 

r 1 

.02  0.519 

0.318 

0.1909 

.04  0.220 

0.1349 

0.0S08 

.06  0.0315 

0.01997 

0.01209 

.08  -0.0815 

-0.0490 

-0.0291 

.10  -0.1424 

-0.0865 

-0.0514 

.12  -0.1687 

-0.1025 

-0.0609 

.14  -0.1718 

-0.1039 

-0.0618 

.16  -0.1599 

-0.0966 

-0.0573 

.18  -0.1381 

-0.0834 

-0.0493 

.20  -0.1107 

-0.0667 

-0.0392 

.30  0,0387 

0.0245 

0.01591 

.40  0.1432 

0.0888 

0.0551 

.50  0.1973 

0.1220 

0.0758 

.60  0.217. 

0.1346 

0.0840 

.70  0.217 

0.1333 

0.0351 

.80  0.210 

0.1316 

0.0834 

1.0  0.1846 

0.1166 

0.0750 

1.2  0.1615 

0.1032 

0.0675 

1.4  0.1434 

0.0926 

0.0614 

1.6 

0.0842 

0.0564 

1.8 

0.0522 

2.0 

2.4 

2.0 

2.4 

2,8 

0.1133 

0.0387 

0.01273 

0.0481 

0.01639 

5.39xlO~3 

7.26xlO"3 

2.44xl0~3 

7.65xlO-4 

-0.01718 

-5.91xlO"3 

-1.99xlO"3 

-0.0304 

-0.01041 

-3.46x10~3 

-0.0360 

-0.01228 

-4.05xlO"3 

-0.0365 

-0.01236 

-4.03xl0"3 

-0.0337 

-0.01133 

-3.64xlO~3 

-0.0290 

-9.58xl0"3 

-2.98x10~3 

-0.0228 

-7.36x10~3 

-2.08xl0~3 

0.01057 

4.95x10"° 

2.51xl0"3 

0.0347 

0.01422 

6.24xlO“3 

0.0477 

0.01965 

8.68xlO~3 

0.0533 

0.0224 

0.01008 

0.0544 

0.0247 

0.01079 

0.0539 

0.0236 

0.01120 

0.0496 

0.0227 

0.01124 

0.0453 

0.0215 

0.01101 

0.0418 

0.0204 

0.01071 

0.0387 

0.01933 

0.01030 

0.0362 

0.01830 

9.86xlO~3 

0.0339 

0.01726 

9 . 40xl0~3 

0.01525 

8.32xl0“3 

7.23xl0“3 

2.8 
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Table  VIII  (Cont.)  Values  of  A'/lfn  for  Cu 


r 

3.2 

3.6 

4.0 

5.0 

6.0 

7.0 

• ^ 
o - 

w 

4.24xl0-3 

1.05xl0-3 

-5.94xlO“6 

1.98xl0“8 

5.93xlO“9 

1.33x10-9 

.04  1.78xl0“3 

4.29xl0"4 

-1.71xlO~3 

2.49xl0“7 

7.45xlO“8 

2.30xl0“8 

.06 

2 . 41xl0-4 

3.3GxlO“5 

-2.69x10“" 

1.07xl0"6 

3.19xlO~7 

9.36xlO“S 

to 

0 

• 

-6.75xl0“4 

-2 . OQxlO“4 

-3.08x3.0~5 

2.86xl0“6 

8.54xlO“7 

2.64xlO“7 

O 

rH 

• 

-I.l6xl0~3 

-3 . 17x10“^ 

-2. 68x10“ 5 

5.93X10"6 

1.78xl0~6 

5.48x10”7 

.12 

-1.33xlO-3 

-3 . 43xl0“4 

-1.38xlO~3 

1.05xl0“5 

3.15xlO~6 

9.73xlO”7 

rH 

• 

-1.31xlO“3 

-3.19xlO“4 

8.47x10“6 

1.67xlO~5 

5.18xlO“6 

1.55xlO“6 

.16 

-1.14xlC“3 

-2.47xl0“4 

3.89x10“5 

2.47xlO~5 

7.37xlO“6 

2.28xl0"6 

to 

rH 

• 

-8.74xl0-4 

-1.45xl0“4 

7.37xlO“5 

3.42xl0”5 

1.03xl0“5 

3.16x10“6 

.20 

-5.51xlO“4 

-1.94x10“"’ 

1. 2^x10“^ 

4.53xlO“5 

1.36xlO~5 

4.18xl0"6 

.30  1.34xlO“3 

7.43x10“^ 

4.30xl0“4 

1.19xl0“4 

3.57xlO~5 

l.lQxlO-5 

.40  2.97xlO"3 

1.49xlO“3 

7.79xlO“4 

2.07xl0“4 

6. 22x10“ 5 

1.92x10“ 5 

.50  4.14xl0“3 

2 . 09xl0~3 

1.10xl0“3 

2.91xl0“4 

8. 74x10“ :> 

2.65xlO“5 

.60  4.91xlO-3 

2.53xlO-3 

1 . 40xlO-3 

3.63xl0~4 

1.09x10“^ 

3.35xlO-3 

.70 

5.39xlO"3 

2.83xl0“3 

1.55xlO“3 

4.18x10“^ 

1.25x10”^ 

3. 87x10“ ^ 

0 

to 

• 

5.69x10”^ 

3.05xl0“3 

1.70xl0“3 

4.62x10"^ 

1.38xl0~4 

4.27xl0-5 

0 

• 

rH 

5.93xl0'3 

3.28xl0"3 

1.36xlO~3 

5 • 18xl0“4 

1.55x10"^ 

4.78x10“5 

1.2 

5.95xlO“3 

3.35xlO“3 

1.94xlO~3 

5.45xlO~4 

1.63xlO~4 

5. 04x10“ 5 

1.4 

5.86xlO“3 

3 .36xlO“3 

1.96xlO~3 

5.54x10-^ 

1.66x10”^ 

5.13xlO“5 

1.6 

5.71xlO“3 

3.29xl0“3 

1.93xlO"3 

5.51xlO“4 

1.65xlO“4 

5. 09x10“ 5 

1.3 

5.51xlO“3 

3.19xlO“3 

1.38xl0“3 

5.39xlO~4 

1.6lxl0“4 

4. 97x10“ ^ 

2.0 

5.27xl0~*3 

3.06xl0“3 

l.SlxlO-3 

5.19xlO~4 

1.56xl0~4 

4. 81x10“ 5 

2.4  4.71xlO“3 

2.74xl0“3 

1.63xlO“3 

4.68x10  4 

1.40xl0“4 

4.33xlO“5 

-3 

2.8  4.11x10 

2.39xl0“3 

1.42xl0~3 

4.10xl0~4 

1.23xl0~4 

3.75xlO“5 
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Table  VIII  (Cont.)  Values  of  A l/Lp  for  Cu 


C\J 

• 

3.6 

4.0 

5.0 

6.0 

7.0 

r 1 

3.2  3.47xl0"3 

2.03xl0"3 

1.21xl0~3 

3.47xl0”4 

l.O/pdO"4 

3.21xlO”3 

3.6 

1.69xl0“3 

1. 00x10 “3 

2.89xl0"4 

8.65xlO“3 

2.67xlO"3 

4.0 

8 . 19xlO”4 

2.36xl0”4 

7. 06x10" 5 

2 . 18x10” 3 

5.0 

1.33xlO“4 

3 . 98x10”^ 

1.23xlO"5 

6.0 

2. 06x10” 3 

6.3AxlO“6 

7.0 

3111xlO"6 

r 8.0 

9.0 

10 

11 

r 1 

.02  5.75xlO"10 

1.94xlO"10 

5.26X10"11 

2.17X10”11 

.04  7.2Axl0”9 

2.44xl0”9 

6.62x10"10 

2.73xl0"lu 

. -8 
.06  3.10x10 

-8 

1.12x10 

2.85xl0-9 

1.17xlO"9 

.08  8.33x10"® 

2.81xlO_S 

7.62xlO"9 

3.13xlO”9 

.10  1.73xlO“7 

5.83x10"® 

1.58x10"® 

6.51xlO"9 

.12  3.06xlO“7 

1.03xl0"7 

2.8QxlO-8 

1.16x10“° 

.14  4.87xlO’7 

1.64xlO~7 

4.46x10 

-8 

1.83x10 

.16  7.17xlO~7 

2.42xl0”7 

6.56xlO-8 

2.70xl0“8 

.18  9.95xl0“7 

-7 

3.36x10 

-8 

9.12x10 

3.75x10”® 

.20  1.32xl0”6 

4.46x10  7 

1.21xl0“7 

4.97x10”® 

* 

-6 

.30  3.47x10 

-6 

1.17x10 

-7 

3.18x10 

1.31xl0“7 

.40  6.02xl0“6 

2 . 04x10”® 

5 . 52xlO-7 

2 . 28xl0~7 

-6 

.50  8.50x10 

-6 

2.86x10 

7.76x10“ 7 

3 . 20xl0“7 

.60  1.15xlO"5 

3.56xlO~6 

9.65x10“ J 

3.98xl0“7 

-5 

.70  1.22x10 

4.11xl0-6 

-6 

1.11x10 

4.58xl0”7 
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Table  VIII  (Concl.)  Values  of  A'/4rr  for  Cu 


r 8.0 9^0 


r'  -5 

.80  1.34x10 

4.51xl0~6 

1.0  1.50xl0~^ 

5.09xl0"6 

1.2  1.58xlO~5 

5.35xlO"6 

-5 

1.4  1.62x10 

5.45xlO~6 

1.6  1.6lxl0“5 

5.42xl0“6 

1.8  1.57xlO"5 

-6 

5.28x10 

2.0  1.51xlO~5 

5 . 10xl0-6 

2.4  1.36x10^ 

4.6QxlO"6 

2.8  1.19x10"° 

-6 

4.02x10 

3.2  l.OlxlO-7 

3.41xlO-6 

, -6 
3.6  8.42x10 

-6 

2.72x10 

, -6 
4.0  6.85x10 

-6 

2.32x10 

5.0  3.S6xl0"6 

-6 

1.30x10 

-6 

6.0  2.00x10 

6.75x10*° 

-7 

7.0  9.79x10 

-7 

3.30x10 

8.0  4.60x10“7 

-7 

1.55x10 

9.0 

-8 

7.49x10 

10 

11 


10  11 


1.23X10""6 

5.08xl0“7 

1.38xlO~6 

-7 

5.69x10 

1.45xlO"6 

5.96xl0"7 

1.48x10-6 

6.09xl0-7 

1.47xlO“6 

6.05xl0“7 

—6 

1.43x10 

5.91xlO-7 

1.38x10 

-7 

5.70x10 

1.25xlO”6 

5.14xlO-7 

-6 

1.09x10 

4.49xl0“7 

9.24xl0“7 

3.81xl0”7 

-7 

7.68x10 

3.17xl0"7 

-7 

6.23x10 

-7 

2.59x10 

-7 

3.53x10 

-7 

1.46x10 

-7 

1.83x10 

-8 

7.54x10 

-8 

8.97x10 

/ ~8 
3.69x10 

-8 

4.21x10 

-8 

1.73x10 

2.02xl0-8 

8.35xlO-9 

7.53xl0“9 

3.10xlO_9 

1.7QxlO-9 
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Table  IX 

Values  of  A* /l^n 

for  Ag 

r .02 

.04 

.06 

.08 

.10 

.12 

r' 

.02  114.8 

88.1 

46.3 

17.23 

0.01225 

-7.83 

.04 

51.0 

37.3 

25.9 

17.58 

11.96 

.06 

37.1 

31.6 

25.7 

20.4 

.08 

32.9 

29.7 

24.3 

.10 

29.8 

24.8 

.12 

23.5 

r .14 

.16 

.18 

.20 

.30 

.40 

r 1 

.02  -11.26 

-10.88 

-9.28 

-7.16 

1.202 

3.17 

.04  7.86 

5.26 

3.48 

2.16 

-0.259 

-0.380 

.06  15.48 

11.43 

8.19 

5.68 

-1.121 

-1.925 

.08  19.09 

14.48 

10.65 

7.36 

-0.899 

-2.12 

.10  19.81 

15.26 

11.47 

8.19 

-0.1774 

-1.655 

.12  19.01 

14.96 

11.54 

8.57 

0.813 

-0.838 

.14  17.74 

14.30 

11.39 

8.88 

1.941 

-0.1466 

.16 

13.52 

11.23 

9.12 

3.09 

1.046 

.18 

11.24 

9.51 

4.18 

1.909 

.20 

9.90 

5.27 

2.76 

.30 

7.80 

4.61 

4.46 
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Table  IX  (Cont.)  Values  of  k'/Ur>  for  Ag 


r .50 

.60 

.70 

r’ 

.02  2.35 

1.095 

0.1458 

.04  -0.213 

-0.0882 

-0.01997 

.06  -1.298 

-0.600 

-0.1048 

.08  -1.492 

-0.695 

-0.1248 

.10  -1.225 

-0.588 

-0.1115 

.12  -0.755 

-0.388 

-0.0845 

.14  -0.1989 

-0.1461 

-0.0433 

.16  0.352 

0.0970 

0.00513 

.18  0.854 

0.321 

0.0552 

.20  1.321 

0.539 

0.1189 

.30  2.43 

1.098 

0.330 

.40  2.43 

1.221 

0.555 

.50  2.03 

1.177 

0.697 

.60 

1.095 

0.801 

.70 

0.867 

o 

to 

• 

1.0 

.80  1.0  1.2 


-0.500 

-0.903 

-0.802 

-3 

4.38x10 

0.01293 

7.76xlO-3 

0.1870 

0.369 

0.324 

0.223 

0.442 

0.-389 

0.1877 

0.372 

0.329 

0.1164 

0.243 

0.218 

0.0424 

0.1027 

0.0947 

-0.0200 

-0.0244 

-0.01750 

-O.C677 

-0.1309 

-0.1112 

-0.0915 

-0.209 

-0.1385 

-0.0395 

-0.278 

-0.259 

0.1925 

-0.0610 

-0.0950 

0.429 

0.1885 

0.1014 

0.617 

0.393 

0.260 

0.744 

0.532 

0.368 

0.811 

0.628 

0.428 

0.658 

0.465 

0.446 

1.2 


I 





Table  IX  (Cont.)  Values  of  A'/Lfn  for  Ag 


r 1.4 

1.6 

1.8 

2.0 

2.4 

2.8 

r' 

.02  -0.582 

-0.405 

-0.264 

-0.1671 

-0.0654 

-0.0257 

.04  4.53xl0"3 

-3 

2 . 66x10 

1.34xlO~3 

8.11xl0*4 

3.25xl0“4 

-4 

1.29x10 

.06  0.234 

0.1623 

0.1033 

0.0666 

0.0259 

0.01018 

.08  0.281 

0.1949 

0.1265 

0.0800 

0.0311 

0.01222 

.10  0.238 

0.1655 

0.1072 

0.0679 

0.0298 

0.01028 

.12  0.1578 

0.1095 

0.0710 

0.0449 

0.01726 

-3 

6.67x10 

.14  0.0692 

0.0404 

0.0302 

0.01949 

7.26xlO"J 

2.67xl0-3 

.16  -0.01153 

-3 

-7.71x10 

-5.22xlO-3 

-3.69xl0“3 

-1.8SxlO"3 

-l.OlxlO-3 

.18  -0.0795 

-0.0547 

-0.0346 

-0.0231 

-3 

-9.55x10 

-4.10xlO~3 

.20  -0.1318 

-0.0909 

-0.0595 

-0.0381 

-0.01549 

-3 

-6.50x10 

.30  -0.1885 

-0.1318 

-0.0842 

-0.0551 

-0.0224 

-3 

-9.32x10 

.40  -0.0757 

-0.0555 

-0.0363 

-0.0234 

-9.53x10^ 

-4.01xlO~3 

.50  0.0611 

0.0374 

0.0244 

0.01578 

-3 

6 . 68x10 

-3 

2.97x10 

.60  0.1854 

0.1128 

0.0738 

0.0481 

0.0205 

9.18xlO~3 

.70  0.247 

0.1639 

0.1078 

0.0706 

0.0307 

0.01396 

.80  0.290 

0.1941 

0.1281 

0.0845 

0.0373 

0.01734 

1.0  0.302 

0.213 

0.1417 

0.0948 

0.0433 

0.0209 

1.2  0.304 

0.206 

0.1385 

0.0938 

0.0442 

0.0220 

1.4  0.280 

0.1909 

0.1294 

0.0887 

0.0429 

0.0219 

1.6 

0.1758 

0.1199 

0.0829 

0.0409 

0.0212 

1.8 

0.1098 

0.0766 

0.0384 

0.0203 

2.0 

0.0707 

0.0360 

0. 01925 

2.4 

0.033-2 

0.01685 

2.8 


0.01439 
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Table  IX  (Cont.)  Values  of  k'/lftr  for  Ag 


r 3.2 

3.6 

4.0 

5.0 

6.0 

7.0 

r ' 

.02  -9.94xl0“3 

-3.95xlO~3 

3.12xl0*"3 

-8.59xl0”e 

-2.59x10”® 

-8.88x10-9 

.04  9.01xl0”5 

2. 02x10”  ■ 5 

7.33xlO”3 

-1.00xl0“6 

-3.02xl0“7 

-9.15x10™® 

.06  3 . 96x10”-'’ 

1.55xlC™3 

9.32xlO“4 

-3.83xlO“6 

-1.15xl0"6 

-3.51xlO“7 

.08  4.73xlO“3 

1.35xl0”3 

7.27xlO_P 

-9.16x10“6 

-2.75xlO"6 

-3.37xlO”7 

.10  3.95xlO~3 

1.59xlO“3 

7.97xlO“5 

-6 

-1.71x10 

-5.14xlO“7 

-1.56xlO“7 

.12  2.50xl0"3 

9.39xl0”4 

-4.54x10“" 

-2.72x10"" 

-8.18xl0-6 

-2.48x10 

.14  3.99xl0“4 

2 . 84xl0”4 

-1.23xl0-4 

-3. 87x10”"’ 

-I.l6xl0~3 

-3.53xlO“6 

.16  -5.7SxlO“4 

-3.20xl0“4 

-2.13xl0~4 

-5. 09x10“ 3 

-1.53x10“" 

-4.65xlO“6 

.18  -1.33xl0~3 

-S.35xl0“4 

-4 

-2.94x10 

-6,26x10”" 

-1. 89x10*""* 

-5.72xlO“6 

.20  -2.81xl0~3 

-1.2ZpclO~3- 

-4 

-3.68x10 

-7. 37x10“ 3 

-9.79X10-6 

-6.73xlO*"6 

.30  -4.00xl0“3 

-1.75xlO-3 

“4 

-5.04x10 

-9.6lxl0“5 

-2.89xlO“5 

—6 

-8.77x10 

.40  -1.75xlO“3 

-7.74xlO-4 

-2.48xl0~4 

-4.73x10"""* 

-1.42x10" 3 

-4.32X10™6 

.50  1.38xl0“3 

6.58xl0“4 

2.69xl0”4 

6.21xl0"3 

1.87xlO“5 

5.68xl0~6 

.60  4.32xl0”3 

2 . 08xl0”3 

8.97xl0“4 

2 . 06xl0“4 

6. 19x10” 3 

1.88xl0“5 

.70  6.73xlO“3 

3 ,32xl0"3 

1.52xl0”3 

3.59xlO“4 

1.08xl0”4 

3.28xlO“3 

.80  8.54xl0~3 

4.30xlO”3 

-3 

2.08x10 

5 . 04xl0~4 

1.52xl0"4 

-5 

4.60x10 

1.0  0.01060 

-3 

5.61x10 

2.91xl0“3 

-4 

7.36x10 

2.21xl0“4 

6.72xl0“5 

1.2  0.01159 

6.25x10“3 

2.55xlO-3 

8.83xl0-4 

2.66x10  4 

S.07xl0-5 

1.4  0.01177 

6 . 48xlO-3 

3.62xlO“3 

9.63xlO“4 

2.89xl0”4 

8.81xl0”3 

1.6  0.01162 

6.49x10~3 

3.69x10“3 

9.96x10-4 

2.99xl0”4 

9.08xl0’3 

1.8  0.01121 

6.33xl0”3 

3.68x10  3 

9.94x10“4 

2.99xl0”4 

9.07xlO"3 

2.0  0.01065 

6.09xl0”3 

3 . 54xlO*"3 

9.71xl0~4 

2.92xl0“4 

-5 

8.87x10 

2.4  9.47xl0“3 

5.43xlO"3 

3.18xl0“3 

8.85xl0~4 

2.66xl0”4 

8.09xl0”5 

2.8  8.15xlO“3 

4.70xl0“3 

2.77xl0”3 

7.72xl0“4 

2.32xl0”4 

7.06xl0“5 
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Table  IX  (Cont.) 


Values  of  A* /Lp  for  Ag 


r 3.2 

r • 

3.2  6.86xlO“3 
3.6 


4.0 


5.0 

6.0 


7.0 


3.6 

4.0 

5.0 

6.0 

7.0 

3.9SxlO~3 

2.35xlO“3 

6.56x10“^ 

1.97xlO“4 

5.99x10-- 

3.28xl0“3 

1.93xl0“'/ 

5.42xl0“4 

1.63xlO"4 

4. 96x10“ 5 

1.57xlO"3 

4.39xl0“4 

1.32xl0”4 

4.06x10~5 

2.40xl0-4 

7.24xl0“5 

2. 20x10“ 3 

3. 75xlO~ 5 1.14x10' 

3.52x10' 


8.0 

9.C 

10 

. i 

02  -2.5QxlO-9 

-7.96xlO”1C 

-3.22xl0“10 

04  -2.93xl0~8 

-9.23xl0“9 

-3.24xl0"9 

06  -1.17x10“ 7 

-8 

-3.54x10 

“8 

-1.43x10 

-7 

08  -2.67x10 

-8.44x10“® 

-3.42x10“® 

10  -4.97xl0“8 

-1.57x10“® 

-6.36xl0“9 

12  -7.96xlO“7 

-2.50xl0“7 

-l.OlxlO-7 

14  -1.13xlO"6 

-3 . 56xlO“7 

-1.44xl0“7 

—6 

16  -1.48x10 

-4.55x10"7 

-1.39xl0"7 

18  -1.84xl0”6 

-5.77xlO"7 

-2.34xl0“7 

20  -2.15x10“® 

-6.79x10“7 

-2.74xl0“7 

-6 

30  -2.80x10 

-3.35xl0~7 

-3 . 59xlO“7 

40  -1.38x10”® 

-4.35xlO“7 

-1.76x10“7 

50  1.81xl0“6 

5.73xlO“7 

2.32xl0“7 

60  5.99xlO“6 

1.89x10“® 

7.68x10“7 

70  1.05xl0-3 

3.31xl0~6 

1.34x10”® 

vn  \_n  vn 
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Table  IX  (Concl.)  Values  o.f  A'/lpr  for  Ag 


r 8.0 

r 1 

.80  1.47x10“ 5 

1.0  2.14x10“" 

1.2  2.57x10”" 

1.4  2.81xl0“5 

1.6  2. 90x10“ ^ 

1.8  2.90x10“'" 

2.0  2.83xl0~5 

2.4  2. 53x10” 5 

2.8  2.25x10“" 

3.2  1.90xlO“5 

3.6  1.5SxlO“5 

4.0  1.28x10“ 5 

5.0  7.02xl0“6 

6.0  3. 64x10  ^ 

7.0  1.76xic"6 

8.0  3.35xlO~7 

9.0 


9.0 

10 

4.65xl0“6 

1.88xl0”6 

6.77xl0~6 

2.74xlO“6 

8.13xlO“6 

3.29xlO”6 

8.85xl0“6 

3.59xlO”6 

9.16x10"6 

3.71x10“° 

9.15xl0-6 

3.71xl0“6 

8.94xl0“6 

3.62xl0“6 

8.14xl0“b 

3.30xl0”6 

7.10xl0"6 

2.88xl0“6 

, -6 
6.04x10 

2.44xl0“6 

5.00xl0“6 

2.02xl0”6 

4.05x10  ~ 

1.64xl0“6 

2.21xl0“6 

8.93x10“7 

1.15xlO“6 

4.66x10~7 

5.55xl0“7 

2.25x10“7 

2.64xlO”7 

1.07xl0“7 

1.19xl0“7 

4.31xlO"S 

, -8 
2.67x10 

10 
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This  process  was  repeated  for  each  r value  of  the  exchange 
integral  needed  in  the  Fock  equation.  The  results  of  these  integrations 
are  shown  in  graphical  form  in  Fig.  6,  where  the  value  of  the  exchange 
integral  is  plotted  as  a function  of  r.  A smooth  curve  was  drawn,  and 
values  actually  used  in  the  solution  of  the  equation  were  taken  from 
this  curve. 

In  the  evaluation  of  the  exchange  integral  for  both  Cu  and  Ag 
it  was  found  that  for  large  values  of  6 the  curve  did  not  follow  the 
trend  of  those  obtained  for  low  6 values.  This  effect  was  not 
pointed  out  by  Fuchs,  nor  by  Gunter,  v.ho  evidently  arbitrarily 
extrapolated  all  curves  to  zero  at  six  atomic  units  distance.  The 
discrepancy  existing  is  readily  noted  by  comparison  of  the  dotted 
line  curves  of  Fig  7 for  an  £ value  of  3.0  with  the  solid  curve 
for  6 = 1.6.  Even  for  € = 1.6  the  curve  tends  to  show  an  increase 
beyond  six  atomic  units,  as  indicated  by  the  dotted  portion. 

It  should  be  remembered  that  the  values  of  A , plotted  as 

r,r«  * * 

ordinates  for  these  curves,  are  obtained  from  the  product  RCr^A* 

r,r» 

Now  the  R- function  values  increase  without  limit  with  increasing  r' . 
However,  the  A’r  r,  values,  as  can  be  seen  from  the  tables,  decrease 
with  increasing  r',  and  are  actually  zero  beyond  eleven  atomic  units. 

Now  for  small  6 values  the  decrease  in  A‘  is  more  rapid 

r,  r' 

than  the  increase  in  R(r’).  Consequently,  the  curves  of  A approach 

r,r« 

the  r1  axis  at  six  or  seven  atomic  units.  However,  for  large  6 values, 
the  increase  in  R(r')  becomes  much  more  rapid  than  the  decrease  in 

A'  , , and  the  A ( values  would  consequently  become  very  great 

a , r r, r 


out  to  eleven  atomic  units,  and  then  suddenly  drop  to  zero. 

Since  the  value  of  the  exchange  integral  at  a given  r is 

represented  by  the  area  of  this  curve,  it  would  not  be  expected  that 

a given  vertical  section  of  this  area,  |Arl,  would  be  greater  at 

r,  r 

nine  atomic  units  than  at  five  atomic  units,  since  the  exchange  forces 
should  be  relatively  small  at  the  former  distance.  It  was  therefore 
concluded  that  an  obvious  error  existed  in  the  functions  as  obtained 
by  Hartree's  method  of  the  self-consistent  field,  and  that  the  true 
shape  of  the  curves  was  approximately  that  given  by  the  use  of  the 
smaller  €■  values.  Consequently,  to  correct  for  this  error,  the  curves 
were  drawn  as  shown  by  the  heavy  lines  in  Fig.  7* 

If  the  areas  of  these  curves,  which  represent  the  value  of  the 
exchange  integral  for  various  r-values,  are  plotted  as  functions  of 
r we  obtain  the  curve  mentioned  previously,  and  shown  in  Fig.  6,  for 
€ = 1.6.  Since  all  of  these  curves  took  the  same  form,  and  approached 
the  r-axis  for  large  r values,  it  was  felt  that  the  correction  for  the 
abnormal  shape  of  the  exchange  integral  curves  was  justified. 


Fig.  7 A2  q r,  vs.  r*  for  Ag  with 
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Figures  8 and  9 show  the  R-function  curves  obtained  without 

exchange  for  copper  and  silver,  respectively.  The  stable  energy 

points  as  given  by  the  condition  dR/dr  = R/r[  are  shown  to  be 

ir  rQ 

nearer  the  nucleus  for  copper  than  for  silver,  as  was  expected. 
Likewise,  a similar  result  is  obtained  with  exchange,  as  shown  in 
Figures  10  and  11.  Figures  12  and  13  show  the  final  E vs.  r curves 
for  both  copper  and  silver.  The  upper  curve,  labled  III,  is  the  final 
curve  obtained  after  adding  the  Fermi  energy  to  the  energy  obtained 
in  curve  II  when  the  differential  equation  was  solved  with  the  exchange 
integral.  Curve  I is  that  obtained  without  exchange.  The  horizontal, 
dotted  line  represents  the  ionization  potential. 

From  these  curves  values  of  the  lattice  energy,  lattice  constant 
and  compressibility  may  be  obtained,  and  compared  with  experimental 
and  previous  theoretical  data  as  shown  in  Table  X. 


Table  X Final  Data  for  Copper  and  Silver 


Theoretical  (as 

Copper 

Cohesive  Energy  Lattice  Constant 
112  Kg.  Cal. /mole  1.98  A 

Compressibility 
0.0082xl0“ij-  cm2 

determined  in  this 
work) 

Experimental 

71.0 

3.6 

dyne 

0.070xl0-11 

Fuchs  (Reported) 

34.1 

4.2 

0.069x10"^ 

(From  his  graphs)  34.1 

4*6 

0.545X10"11 

Theoretical  (as 

Silver 

71.8  Kg.  Cal. /mole 

2.57  A 

0.035xl0“1]-  cm2 

determined  in  this 
work) 

Experimental 

64.5 

4.03 

dyne 

0.090X10"11 

Gunter 


95.8 


2.72 


(0.056x10_11) 


1 


f 


Atomic  Units 


45 


46 


r.  Atomic  Units 
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Rydberg  Units 
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Conclusions 


An  examination  of  the  final  data  presented  in  Table  X brings 
out  several  interesting  points.  In  the  case  of  silver,  agreement 
between  theoretical  and  experimental  values  for  this  element  are  much 
better  than  for  copper.  The  writer  obtains  a slightly  lower  value 
for  the  lattice  constant  of  silver  than  does  Gunter,  his  value 
being  closer  to  the  experimental  value,  but,  on  the  other  hand,  a 
much  better  value  for  the  cohesive  energy  is  obtained  in  this 
investigation.  The  reasons  for  the  differences  that  do  exist  are 
not  self  explanatory,  but  probably  are  associated  with  the  methods 
used  in  solving  the  differential  equations  and  the  extent  of  the 
radial  wave  functions  used  in  evaluating  the  exchange  integral. 

Since  Gunter  did  not  record  a value  for  the  compressibility,  the 
value  he  would  have  obtained  from  his  energy  curve  is  shown  in 
parentheses. 

In  the  case  of  copper,  it  is  difficult  to  tell  wherein  the  error 

lies.  The  discrepancy  existing  between  values  reported  by  Fuchs  and 

those  obtained  by  calculation  from  his  graph  is  not  understandable. 

On  the  other  hand,  the  close  agreement  the  writer  has  with  the  work 

of  Gunter  on  silver,  and  the  fact  that  the  changes  in  values  as  one 

goes  from  silver  to  copper  are  in  the  proper  direction,  lead  the  writer 

to  feel  that  his  results  are  the  more  acceptable.  Since  Fuchs'  value 

for  the  compressibility  is  actually  0.545  cm  /dyne  rather  than  his 

2 

reported  value  of  O.O69  cm  /dyne,  there  is  no  reason  to  place  emphasis 
on  values  obtained  for  this  physical  property. 
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The  ion-ion  repulsions  were  not  added  in  this  work,  since  the 
amount  of  labor  involved  in  their  calculation  in  addition  to  what 
has  already  been  done  would  be  prohibitive.  The  addition  of  the 
ion-ion  repulsion  should  give  better  results  for  all  quantities, 
however,  providing  it  would  not  be  too  large.  Any  values  for  this 
repulsion  energy  near  those  reported  by  Fuchs  would  be  much  too  great 
to  apply  as  corrections  here,  and  it  is  possible  that  a completely 
new  theoretical  treatment  of  ion-ion  repulsion  is  in  order. 

The  results  obtained  thus  far  indicate  an  error  of  major 
importance  in  the  previous  calculations  for  copper.  It  is  hoped 
that  this  work  vail  open  anew  the  subject  of  cohesive  energies  and 
lattice  constants  for  the  noble  metals  with  a view  to  obtaining  a 
complete  set  of  consistent,  theoretical  data  for  copper,  silver  and 
gold. 
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Comprehensive  Abstract 


E Wigner  and  F.  Seitz  presented,  perhaps,  the  first,  simplified 
quantum  mechanical  treatment  of  cohesive  energy,  lattice  constant 
and  compressibility  of  metallic  crystals  in  their  work  on  the 
constitution  of  metallic  sodium.  They  applied  the  free  electron 
assumption,  considering  that  the  inner,  core  electrons  were  not 
affected  by  the  metallic  bond,  and  that  a satisfactory  treatment 
could  be  given  by  applying  the  radial  form  of  the  Schroedinger 
equation  to  the  valence  electrons.  The  boundary  condition  was 
obtained  by  assuming  that  the  V^)  -function  would  be  a repeating 
function  in  going  from  atom  to  atom  in  the  lattice,  and  that  the 
derivative  of  the  function  at  every  crystallographic  symmetry  plane 
would  be  zero  perpendicular  to  this  plane.  Thus  if  the  -function 
is  such  that  R = r Vp  , where  R is  the  radial  function  and  r is  the 


distance  out  from  the  center  of  the  nucleus  the  boundary  condition 

<A*Va 

^ so  that  dR/dr  = R/r. 


is  given  by 


- o 


This  condition  would  occur  on  the  R-curve  at  points  where  a 
tangent  would  cut  through  the  origin.  The  point  of  tangency  would 
then  give  the  r value  for  the  most  stable  energy  point.  In  general, 
every  R-curve  for  energy  less  (numerically)  than  that  for  the  most 
stable  energy  position  would  give  two  points  of  tangency.  For  energy 
greater  (numerically)  there  would  be  no  point  of  tangency. 

Using  this  procedure,  applying  potential  functions  worked  out 
by  Prokof jew,  and  considering  the  Fermi  energy  and  correlation  energy, 
Wigner  and  Seitz  obtained  a value  of  23.2  kilogram  calories/mole  as 


c 
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against  the  experimental  value  of  26.9  kilogram  calories/mole  for  the 
cohesive  energy.  For  the  lattice  constant  they  obtained  4.75  A as 
compared  with  4.23  A obtained  experimentally. 

A similar  method  was  applied  to  lithium  by  Seitz,  who  obtained 
33.3  kilogram  calories/mole  as  against  38*9  kilogram  calories/mole 
for  the  cohesive  energy,  and  3.53  A as  compared  with  the  experimental 
value  of  3.46  A for  the  lattice  constant. 

E.  Gorin  attempted  a theoretical  consideration  of  potassium  by 
this  method,  but  found  that  he  could  not  employ  Prokof jew  potentials 
with  any  degree  of  success.  Consequently,  he  eventually  used  a set 
of  Hartree  functions  for  all  of  the  electrons,  and  evaluated  the 
exchange  energy  directly.  He  obtained  14.3  kilogram  calories/mole 
as  against  22.6  kilogram  calories/mole  for  the  experimental  value  of 
the  cohesive  energy. 

The  first  attempt  at  the  application  of  Hartree  functions  to 

the  solution  of  a Fock  equation  for  one  of  the  noble  metals  was  made 

by  K.  Fuchs.  He  applied  a set  of  available  Hartree  functions  for 

copper,  and  evaluated  the  exchange  integral,  directly.  After  adding 

the  Fermi  energy  he  reported  values  of  41.4  kilogram  calories/mole, 

-11  2 

4.1  A,  and  0.26  x 10  cm^/dyne  for  the  cohesive  energy,  lattice 

constant,  and  compressibility^,  respectively.  These  might  be  compared 

with  corresponding  experimental  values  of  71.0  kilogram  calories/mole, 

-11  2 

3.6  A,  and  0.07  x 10  cm  /dyne. 

Fuchs  then  claimed  to  obtain  a better  value  for  the  compressibility 
by  considering  the  ion-ion  repulsion.  VJhen  this  energy  is  added  he 


* 


i 


K. 


t 


C 


—11  2 

reports  34.1  kilogram  calories/mole,  4.2  A,  and  0.069  x 10  cm  /dyne, 
respectively,  for  these  values.  A criticism  of  these  results  together 
with  a completely  new  set  of  calculations  forms  a major  portion  of 
the  dissertation. 

It  appears  that  nothing  more  was  done  in  this  field  following 
Fuchs'  work  until  R.  C.  Gunter  attacked  a similar  problem  for  silver. 
Ke  applied  exactly  the  same  method  that  Fuchs  did  for  copper,  a set 
of  Hartree  functions  having  been  available.  The  Fermi  energy  was 
added,  and  he  reported  values  of  95  • 8 kilogram  calories/mole  and 
2.72  A for  the  cohesive  energy  and  lattice  constant,  respectively, 
the  experimental  values  being  64.5  kilogram  calories/mole  and  4.03  A. 

A careful  study  of  the  energy  curves  obtained  for  copper  by- 
Fuchs  and  for  silver  by  Gunter  shows  the  stable  energy  point  to  be 
higher  for  copper  than  for  silver.  Obviously,  if  the  same  method  had 
been  used  for  both  atoms,  one  of  these  results  must  be  in  error. 

Also,  Fuchs  reports  that  the  reason  for  considering  the  ion-ion 
repulsion  was  to  obtain  a better  value  for  the  compressibility. 
Actually,  the  value  he  would  have  obtained  without  considering  the 
ion-icn  repulsion  energy  gives  a larger  value  for  r^,  the  stable 
■energy  point,  and  hence  a higher  value  for  the  compressibility. 

In  view  of  this  discrepancy,  the  writer  decided  to  work  through 
independently  the  wave  functions  without  and  with  exchange  for  both 
copper  and  silver,  and  determine,  if  possible,  the  source  of  error. 

In  the  course  of  the  investigation  it  was  necessary  to  solve  two 
differential  equations  and  one  integral  by  numerical  methods  for  both 
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copper  and  silver.  The  equations  to  be  solved  were  forms  of  the 


Schroedinger  equation,  the  one  without  exchange  being  simply 

♦ ^(E  - V)B  - 0 

where  R = r^(r),  r being  the  distance  out  from  the  center  of  the 
atom,  and  y being  the  eigenfunction  for  the  valence  electron. 
Using  atomic  units  for  length  and  energy: 

1 atomic  unit  of  length  = h2/Arr2me2  = 0.528  A 
1 atomic  unit  of  energy  = 2tt  meH  = 1 Rydberg  unit 
1 Rydberg  unit  = 13.53  electron  volts/atom 
= 312  kilogram  calories/mole 
If  we  also  set  E = - C/2,  the  equation  becomes 

d2l 

-JP?  + (2V  - € )R  = C 


which  is  the  form  in  which  it  is  used  for  determining  the  R-function 
without  exchange. 

The  equation  to  be  solved  with  exchange  takes  the  form 


d2R 


^ + (2V  - 6)R  + 2 [G(r,r')R(r,)dr'  = 0 
dr  J0 

I G(r  r 1 )R(r’ )dr 1 is  the  exchange  integral. 

Jo 

into  the  form  e2  ^ K0f  w Ir)fw  „(r* ) , where 

K 0 * K.-T  ’ 


where 
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G(r, r ' ) may  be  put 
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The  f's  are  the  Hartree  functions  for  the  individual  electrons; 
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thus  0(0.02)  would  represent  the  value  of  the  radial  wave  function 
for  a 3s  electron  at  a distance  of  0.02  atomic  unit  from  the  center  of 
the  atom.  R(r')  is  the  value  of  the  R-function  to  be  inserted  from 
the  solution  of  the  equation  without  exchange.  The  value  of  the  integral 
for  a given  r is  inserted  in  the  equation  together  with  a corresponding 
value  of  V as  obtained  from  the  Hartree  field.  6 is  approximated 
until  a value  is  found  for  which  there  is  only  one  tangent  to  the 
R-curve  as  obtained  with  exchange.  This  will  represent  the  r-position 
of  the  most  stable  energy  point.  The  Fermi  energy,  or  mean  kinetic 
energy,  is  then  added.  This  is  obtained  from  the  following  expression: 

-5  f 3 1 _ 'g--'z-gm 

vo \l^  ~ n* 

3 

V is  the  atomic  volume.  4/3  trr  • 

If  EQ  is  the  self  energy,  that  is,  the  energy  without  exchange, 
of  the  valence  electron,  E„  the  Fermi  energy,  and  E^  the  negative  of 
the  ionization  energy,  we  can  express  the  cohesive  energy  as 

Ec  = E0  + Ef  + Ei  • 

The  cohesive  energy  is,  then,  the  difference  between  the  ionization 
energy  and  the  energy  the  electron  has,  including  Fermi  energy,  when 
the  atom  is  in  a lattice  structure. 

The  compressibility  was  calculated  from  the  formula 

K ■ ^(W3)l/V 


It  should  be  pointed  out  that  experimentally  determined  values  of 
lattice  constant  and  compressibility  would  not  follow  this  curve, 


. 


c 


5 


c < 


t 


< 


i 


hence  although  a lattice  constant  value  determined  theoretically 
would  correspond  with  the  experimentally  determined  value,  the 
compressibilities  would  not  in  general  correspond.  Potassium  is  the 
only  one  of  the  alkalis  or  noble  metals  whose  experimental  values  of 
compressibility  and  lattice  constant  fall  on  this  curve.  The  lattice 
constant,  incidentally,  is  related  to  the  most  stable  energy  point, 
rQ,  by  the  expression 

a = (l6n/3)1/3(0.52S)ro 

This  gives  the  value  in  Angstrom  units  when  Tq  is  in  atomic  units, 
the  factor  0.528  being  the  conversion  factor. 

The  results  of  the  investigation  compared  with  corresponding 
experimental  and  previous  theoretical  data  are  shown  in  the  following 
table: 


Copper 


Theoretical  (as 
determined  in  this 
work) 

Experimental 

Cohesive  Energy  Lattice  Constant 
112  Kg.  Cal. /mole  1.98  A 

71.0  3.6 

Compressibility 
0.0082X10-11  cm2 

dyne 

0.070X10'11 

Fuch3  (Reported)  34.1 

(From  his  graphs)  34.1 

4.2 

4.6 

0.069xl0-11 

0.545xl0“H 

Theoretical  (as 
determined  in  this 

Silver 

71.8  Kg.  Cal. /mole 
work) 

2.57  A 

0.035X10”11  cm2 
dyne 

Experimental 

64.5 

4.03 

0.090xl0-11 

Gunter 

95.8 

2.72 

(0.056x10~11) 

An  examination  of  these  data  brings  out  several  interesting  points. 
In  the  case  of  silver,  agreement  between  theoretical  and  experimental 
values  for  this  element  are  much  better  than  for  copper.  The  writer 


t 


obtains  a slightly  lower  value  for  the  lattice  constant  of  silver 
than  does  C-unter,  his  value  being  closer  to  the  experimental  value, 
but,  on  the  other  hand,  a much  better  value  for  the  cohesive  energy 
is  obtained  here,  ihe  reasons  for  the  differences  that  do  exist  are 
not  self  explanatory,  but  probably  are  associated  with  the  methods 
used  in  solving  the  differential  equations  and  the  extent  of  the  radial 
wave  functions  used  in  evaluating  the  exchange  integral. 

In  the  case  of  copper,  it  is  difficult  to  tell  wherein  the  error 

lies.  The  discrepancy  existing  between  values  reported  by  Fuchs  and 

those  obtained  by  calculation  from  his  graph  do  not  correlate.  On  the 

other  hand,  the  close  agreement  the  writer  has  with  the  work  of  Gunter 

on  silver,  and  the  fact  that  the  changes  in  values  as  one  goes  from 

silver  tc  copper  are  in  the  proper  direction,  lead  the  writer  to  feel 

that  his  results  are  the  more  acceptable.  Since  Fuchs'  value  for  the 

compressibility  is  actually  0.545  x 10  cm  /dyne  rather  than 
-11  2 

0.069  x 10  cm  /dyne,  as  he  reports,  there  is  no  reason  to  place 
emphasis  on  values  obtained  for  this  property. 

The  ion-ion  repulsions  were  not  added  in  this  work,  since  the 
amount  of  labor  involved  in  their  calculation  in  addition  to  what  has 
already  been  done  would  be  prohibitive.  The  addition  of  the  ion-ion 
repulsion  should  give  better  results  for  all  quantities,  however, 
providing  it  would  not  be  too  large.  Any  values  for  this  repulsion 
energy  near  those  reported  by  Fuchs  would  be  much  too  great  to  apply 
as  corrections  here,  and  it  is  possible  that  a completely  new 
theoretical  treatment  of  ion-ion  repulsion  is  in  order. 
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The  results  obtained  thus  far  indicate  an  error  of  major 
importance  in  the  previous  calculations  for  copper.  It  is  hoped 
that  this  work  will  open  anew  the  subject  of  cohesive  energies  and 
lattice  constants  for  the  noble  metals  with  a view  to  obtaining  a 
complete  set  of  consistent,  theoretical  data  for  copper,  silver 
end  gold. 
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